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Abstract. This final report summarizes the work funded under the Grant NAG3-2411 during 
the 04/05/2000-04/04/2001 period. The objective of this one-year project was to generalize the 
theoretical framework of the two-dimensional higher-order theory for the analysis of cylindrical 
functionally graded materials/structural components employed in advanced aircraft engines devel- 
oped under past NASA-Glenn funding. The completed generalization significantly broadens the 
theory’s range of applicability through the incorporation of dynamic impact loading capability into 
its framework. Thus it makes possible the assessment of the effect of damage due to fuel impurities, 
or the presence of submicron-level debris, on the life of functionally graded structural components. 
Applications involving advanced turbine blades and structural components for the reusable-launch 
vehicle (RLV) currently under development will benefit from the completed work. The theory’s 
predictive capability is demonstrated through a numerical simulation of a one-dimensional wave 
propagation set up by an impulse load in a layered half-plane. Full benefit of the completed gener- 
alization of the higher-order theory described in this report will be realized upon the development 
of a related computer code. 


1. INTRODUCTION 

Functionally graded materials (FGMs) are a new generation of composites wherein the microstrue- 
tural details are spatially varied through nonuniform distribution of the reinforcement phase(s), by 
using reinforcement with different properties, sizes and shapes, as well as by interchanging the roles 
of reinforcement and matrix phases in a continuous manner. The result is a microstructure that 
produces continuously changing thermal and mechanical properties at the macroscopic or contin- 
uum level. This new concept of engineering the material’s microstructure allows, for the first time, 
to fully integrate both the material and structural considerations into the final design of structural 
components. 

Most computational strategies for the response of FGMs do not explicitly couple the material’s 
heterogeneous microstructure with the structural global analysis. Rather, local effective or macro- 
scopic properties are first obtained through homogenization based on a chosen micromechanics 
scheme, and then used in a global thermomechanical analysis. This often leads to erroneous results 
when the temperature gradient is large with respect to the dimension of the inclusion phase, the 
characteristic dimension of the inclusion phase is large relative to the global dimensions of the com- 
posite, and the number of uniformly or nonuniformly distributed inclusions is relatively small (Ref. 
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[1]). As a result of the limitations of the uncoupled approach, a new higher-order micromechanical 
theory for FGMs (HOTFGM), that explicitly couples the local and global effects, has been devel- 
oped in the Cartesian coordinate system for applications involving rectangular plate-like structural 
components under NASA funding (Refs. [2-19]). The development of the theory has been justified 
by comparison with the results obtained using the standard micromechanics approach which ne- 
glects the micro-macrostructural coupling effects (Refs. [4,7]). Summaries of significant results and 
accomplishments generated using this theory, and the utility of functionally graded microstructures 
in enhancing the performance of plate-like structural components subjected to through-thickness 
thermal gradients have been outlined in Refs. [9,18,19,20]. 

In order to exploit the already-proven predictive capabilities of HOTFGM to the fullest in meet- 
ing the challenges and needs of the aerospace and aircraft engine industries for a greater number 
of stronger, lighter and more durable structural components, new versions of HOTFGM have been 
recently developed for applications involving cylindrical bodies of revolution. The generalization of 
HOTFGM to problems involving cylindrical geometries makes possible the analysis, optimization 
and design of functionally graded structural components, such as rotor disks, combustor linings and 
blisk blades, for usepn advanced aircraft engines. Thus far, two versions of the cylindrical higher- 
order theory have been developed. The quasi one-dimensional version enables analysis, design and 
optimization of cylindrical bodies of revolution subjected to axisymmetric thermomechanical load- 
ing that are reinforced by either continuous or discontinuous fibers with variable spacing in the 
radial direction (Ref. [21]). The recently completed development of the two-dimensional version, 
which also admits the presence of cooling channels, enables analysis, design and optimization of 
cylindrical bodies of revolution with functionally graded microstructures in the radial and circum- 
ferential directions (Ref. [22]). The thermomechanical loading involves arbitrary distribution of 
surface tractions and temperatures applied to the boundaries of fully or partially enclosed cylin- 
drical bodies of revolution in the plane that contains the functionally graded microstructure. This 
significant generalization provides analysis and design capabilities for a wider range of structural 
components employed in advanced aircraft engines. 

The utility of the developed Cartesian and cylindrical versions of the higher-order theory has 
been demonstrated through applications to the following technologically important problems: 

• Investigation of the effect of microstructure on thermal and stress fields in MMC plates and 
cylinders 

• Investigation of the use of functionally graded architectures in reducing edge effects in MMC 
plates 

• Optimization of functionally graded microstructures in MMC plates and cylinders 

• Development of guidelines for the design of special coatings in exhaust nozzle applications 
under NASA/Pratt & Whitney Space Act Agreement 

• Investigation of the microstructural effects in functionally graded TBCs 

• Effect of bond coat interfacial roughness amd oxide film thickness on the inelastic response 
of plasma-sprayed TBCs 

• Effect of graded bond coats on the inelastic response of plasma-sprayed TBCs 
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While the recently completed cylindrical higher-order theory contains transient thermal loading 
capability, the mechanical loading capability does not include dynamic effects appearing in the 
governing force equilibrium differential equations. This, in turn, excludes the possibility of gaging 
the effect of impact loading (by fuel impurities or submicron debris, for instance) on damage evo- 
lution in such applications as advanced, functionally graded thermal barrier coatings, for instance. 
The objective of the work summarized in this report, therefore, was to extend the two-dimensional 
higher-order theory for cylindrical functionally graded structural components by incorporating dy- 
namic impact loading capability. This extension funded under the Grant NAG3-2411 during the 
04/05/2000-04/04/2001 period makes possible the assessment of the effect of foreign object impact 
on the potential for damage evolution in advanced turbine blade coatings. It complements current 
nation-wide efforts by a number of government agencies to develop a new generation of turbine 
blade coatings capable of operating longer in low-cost fuel environments containing different types 
of impurities. As the impact problem also plays an important role in a number of other techno- 
logically important applications that are important to the nation’s security interests (graded body 
armour, for instance), the completed work significantly broadens the range of technologically im- 
portant applications of the cylindrical higher-order theory. However, the completed work is limited 
to the development of the theoretical framework that enables modeling of impact-induced wave 
propagation in the radial and circumferential directions of functionally graded cylindrical struc- 
tural components. Therefore, the description of this theoretical development forms the major part 
of this report. Full utilization of the impact-loading capability requires the development of the 
related computer code and its validation. These tasks remain to be completed under future fund- 
ing. To demonstrate the potential of the developed theoretical framework, a small computer code 
was developed to simulate one-dimensional wave propagation due to impulse loading in a layered 
half-plane as a special case of the general two-dimensional theory, and the numerical predictions 
were validated through comparison with an exact analytical solution. These results are described 
at the end of the report. 


2. MODEL DESCRIPTION 

The present FGM theory is based on the geometric model of a heterogeneous composite occupying 
the region Rq < r < Ru o < 9 < ©, | 2 |< oo, where r,9,z are cylindrical coordinates, see Fig. 1. 
The composite is reinforced in the r — 9 plane by an arbitrary distribution of infinitely long fibers 
oriented along the axial 2-axis, or by finite-length inclusions that are arranged in a periodic manner 
in the axial direction. The microstructure of the heterogeneous composite is discretized into N p 
and N q cells in the intervals Rq <r < R\ and 0 < 9 < ©, respectively. As in the Cartesian version 
of the higher-order theory (see Ref. [19], for instance), the generic cell ( p,q,s ) used to construct 
the composite consists of eight subcells designated by the triplet (a/3 7), where each index a,jd , 7 
takes on the value 1 or 2 to indicate the relative position of the given subcell along the r — , 6—, and 
z— axis, respectively. The indices p and q , whose ranges are p = 1,2, ...,N P and q = 1,2 ,...,N q , 
identify the generic cell in the r — 6 plane and thus remain constant along the axial z axis. For the 
axial direction, the corresponding index s having an infinite range is introduced. The dimensions 
of the generic cell along the periodic axial direction, luh, are fixed for the given configuration; 
whereas the dimensions along the r— and 6— axes or the FG directions, af , af' , and 9f ,9f , can 
vary in an arbitrary fashion such that D = + df^) and © — Yl^=i(9[ q + 9^)- 

Given the applied thermomechanical loading in the ( r — 0)-plane, an approximate solution for 
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the temperature and the time-dependent displacement fields is constructed based on the volumetric 
averaging of the field equations together with the imposition of boundary and continuity conditions 
in the average sense between the subvolumes used to characterize the material’s microstructure. 
This is accomplished by approximating the temperature field in each subcell of the generic cell 
using a quadratic expansion in the local coordinates (f^, 8^\ z^) centered at the subcell’s center. 
Similarly, the time-dependent displacement field in the FG direction in each subcell is approximated 
using a quadratic expansion in local coordinates within the subcell. The displacement field in 
the periodic axial direction, on the other hand, is approximated using linear expansion in local 
coordinates to reflect the periodic character of the composite’s microstructure along the axis. A 
higher order representation of the temperature and displacement fields is necessary to capture the 
local effects created by the thermomechanical field gradients, the microstructure of the composite, 
and the finite dimensions in the FG direction. 

The unknown coefficients associated with each term in the temperature field expansion are then 
obtained by constructing a system of equations that satisfies the requirements that the steady state 
heat equation is satisfied in a volumetric sense, and the thermal and heat flux continuity conditions 
within a given cell, as well as between a given cell and adjacent cells, are imposed in an average 
sense, together with the applied boundary conditions. 

Due to the presence of the inertia effects in the governing mechanical equations, on the other 
hand, the second time derivatives of the unknown time-dependent coefficients associated with each 
term in the displacement field expansion are obtained in this case by constructing a system of equa- 
tions that satisfies the requirements that the elastodynamic equations are satisfied in a volumetric 
sense, and the displacement and traction vectors continuity conditions within a given cell, as well as 
between a given cell and adjacent cells, are imposed in an average sense, together with the applied 
time-dependent boundary conditions. This system of second order ordinary differential equations 
is solved in a stepwise manner in time by employing an approximate explicit scheme. Once these 
coefficients have been determined at the current time all field variables can be readily established. 
This procedure is continued until the specified final time is reached. 

3. THERMAL ANALYSIS 


Let the composite be subjected to the steady-state temperature distributions Tt(6) on the top 
surface (r = R\), Tb(0) on the bottom surface (r — Rq), T^(r) on the left surface (6 = 0), and 
Tft(r) on the right surface (6 = 0). Under these circumstances, the heat flux field in the material 
occupying the subcell (a/?q) of the (p , q , s)th cell must satisfy the steady state heat equation in 
cylindrical coordinates (r,9,z)> This equation is given by 


H 


(ctPi) 


Of {a) 


,(a/37) 


+ 


n (ccj3y) 

()y(P) 


(afty) 


+ 


di 

dzi'r'i 


= 0 


( 1 ) 


where R( a ^ is the distance of the (afty) subcell’s center from the origin, = R^ a ^6^\ and 
^(a/3-y) ( i — r t 9,z) are the components of the heat flux vector in the subcell. These components 
are derived from the temperature in the subcell according to the Fourier law: 


— _k( Q 07) 
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(a/37) _ i 

q ° ~ 8 dyW 


q{ap~t) _ _£.(a/ 37) 


QJ'(a0l) 

dzW 


( 3 ) 

( 4 ) 


where /cj a ^ denote the thermal conductivities of the material in the subcell. 

As in the analysis of the steady-state heat equation in Cartesian coordinates, the temperature 
field in the subcell is expanded quadratically in terms of the local coordinates as 

follows 


rn(afiy) rp(a /3 7) -(a)o-i(a/?7) . -(ff) rp (or/3-y) 

x x (C\C\(\\ T f - i (l00) ' * -‘Ynim 
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1 ^(p) 2 
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J J (200) 


/l 
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(«) 2 1 /2 
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( 5 ) 


where = R^ a ^6^\ The microvariables Tj^j\ which is the volume- averaged temperature 

within the subcell, and (^ m , n = 0, 1, or 2 with l + m + n < 2) are unknown coefficients 

that axe determined in the manner described below. It should be noted that the temperature 
expansion given in eqn (5) does not contain a linear term in the local coordinates z ^ . This follows 
directly from the assumed periodicity in the axial z-direction and symmetry with respect to the 
lines = 0 for 7 = 1 and 2. 

Given the six unknown quantities associated with each subcell (i.e. —1^002^) an< 4 

eight subcells within each generic cell, 48N p N q unknown quantities must oe determined for a 
composite with with N p and N q subcells containing arbitrary specified materials. These quantities 
axe determined by first satisfying the heat conduction equation, as well as the first and second 
moment of this equation in each subcell, in a volumetric sense in view of the above temperature 
field approximation. Subsequently, continuity of heat flux and temperature is imposed in an average 
sense at the interfaces separating adjacent subcells as well as neighboring cells. Fulfillment of these 
field equations and continuity conditions together with the imposed thermal boundary conditions at 
the top, bottom, left and right surfaces of the composite provides the necessary 48 N p N q equations 
for the 48N p N q unknown coefficients in the temperature field expansion. We begin the outline of 
steps to generate the required 48 N p N q equations by first considering an arbitrary (p, q, s)th cell 
in the interior of the composite (i.e., p = 2,...,N P — 1 and q = 2 — 1). This produces 
48(IVp - 2 )(N q - 2) equations. The additional equations are obtained by considering the boundary 
cells (i.e., p = l,jV p and q = 1 ,N q ). For these cells, most of the preceding relations also hold, 
with the exception of some of the interfacial continuity conditions between adjacent cells that are 
replaced by the specified boundary conditions. 

In the course of satisfying the steady-state heat equation in a volumetric sense, it is convenient 
to define the following flux quantities: 


[Q 


(a/37) 



s ) 


— 1 — rdn rh</2 fV* (a)t .W ) - ( ^) ) n>^) df (a)^W df(7 ) 

(p,9,s) J-JM / 2 J-hfl 2 J-ly/ 2 V 

(a/> 7) p 


(6) 


where i = r, 0,z; Z,m, n = 0, 1, or 2 with l + m + n < 2 and being the volume 

of the subcell. For l = m = n = 0, <2-“ 0 ^o) > s the average value of the heat flux component qf^ 
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in the subcell, whereas for other values of ( l,m,n ) this equation defines higher order heat fluxes. 
These flux quantities can be evaluated explicitly in terms of the coefficients by performing 

the required volume integration. This yields the following nonvanishing zeroth- and first-order heat 
fluxes in terms of the unknown coefficients in the temperature field expansion 


'‘Mo, o,o) 
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if 

(10) 

]2 

" ** 4 ± (002) 
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Satisfaction of the zeroth, first, and second moment of the steady-state heat equation (1) results 
into the following eight relationships among the first-order heat fluxes in the different 

(a/? 7 ) subcells of the (p, q, s)th cell, after some involved algebraic manipulations (see Ref. [12] for 
a complete derivation in Cartesian coordinates) 


r r\( a P'f) , 12 n (a/3y) 
l d 2^r(l,0,0)+ 


+ —Q 

0 ( 0 , 1 , 0 ) T ;2 


(a/37) 


2(0,0, 1) + £(<*/ 97) ^r(0,0,0) 


Q%U p ' qj) = 0 


( 12 ) 


where the triplet (a(3y) assumes all permutations of the integers 1 and 2. 

The continuity of the heat fluxes at the subcell interfaces and between individual cells in the 
radial direction, imposed in an average sense, is ensured by 
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The coefficients Zi, i = 1, ...8, in eqns (13)-(14) are defined as follows 
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The continuity of the heat fluxes at the subcell interfaces and between individual cells in the 
0-direction, imposed in an average sense, is ensured by 




■hi 

^(q27) 




hV 

(a27) 


(17) 


- +«[pQ^! 0) ] c 

2 2 

The coefficients Y u i = 1, ... 8 , in eqns (17)-(18) are defined as follows 
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Finally, the continuity of the heat fluxes at the subcell interfaces in the z-direction, imposed in 
an average sense, is ensured by 


r-L/o (a/31) 4. 

1^(0, o,i)+ i 2 


l^(a/32) i(p, 9 ,s) _ q 

[o V*(0 1 0,i)i - u 


( 21 ) 


The above equations (13)-(14), (17)-(18) and (21), provide us with 20 additional relations 
among the zeroth- and first-order heat fluxes. These 28 relations can be expressed in terms of the 
unknown coefficients by making use of the expressions for heat fluxes given in terms of these 

coefficients in eqns (7)-(ll). 

An additional set of 20 equations that are necessary to determine the unknown coefficients in 
the temperature field expansion is subsequently generated by the thermal continuity conditions 
imposed on an average basis at each subcell and cell interface. Imposing the thermal continuity at 
each subcell interface and between individual cells in the r-direction we obtain 
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In the 0-direction we have 
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The thermal continuity conditions in the periodic z-direction, imposed in the average sense, 
provide 


Trp{QLPl) l\ rpjctP 1) rp(a/32) ^2 rp(a02) 1 (p,Q,s) A 

\r ( 000 ) ^ 4 1 ( 002 ) 1 ( 000 ) 4 1 ( 002 ) J ~ U 
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These temperature continuity conditions, eqns (22)-(26), comprise the required additional 20 rela- 
tions. 

The steady-state heat equations (12) together with the heat flux, (13)-(14), (17)-(18), and 
thermal continuity, (22)- (26), equations form altogether 48 linear algebraic equations that govern 
the 48 field variables in the eight subcells (a/3j) of the interior cell (p,g, $),p = 2, ...N p —l,q = 

2, ...)N q — 1. For the boundary cells p = 1 , N p and q = 1 , a different treatment must be applied. 

For p = 1, the flux continuity conditions (13)-(14) between a given cell and the preceding one are 
not applicable. They are replaced by the condition that the heat flux at the interface between 
subcells (l/? 7 ) and ( 2 /?q) of cell (l,g, s) is continuous as well as the applied temperature relation 
at the surface r = Rq. For the cell p = ' N pi the previous equations are applicable except for 
those which express continuity between this cell and the next one, eqn (23), These equations are 
replaced by the boundary conditions that are applied at the surface r = R\. In the case in which 
the temperature is prescribed, the boundary conditions at the bottom and top surfaces are 


t (i/? 7 ) |(l, <?.*)= t b (9) f (1) = -<4 1) /2 (27) 

t ( 2 /? 7 ) | (N p ,q,s) = Tt ^ f ( 2 ) _ j 2 (28) 

where q — 1, N q , 

Similarly, continuity conditions (17)-(18) that are not applicable at 0 ~ 0 at cell (p, 1 , s) are 
replaced by the continuity of heat flux at the interfaces between the subcells of this cell, and by 
the applied loading at the surface 0 = 0. The temperature continuity conditions between a cell and 

the next one in the ^-direction, eqn (25), which are not applicable at cell (p, N q , s) are replaced by 

the applied loading conditions at 0 = 0. In the case in which the temperature is prescribed, the 
boundary conditions at the left and right surfaces are 

r (ai7) x L (r) y w = -h^/2 (29) 

T (a2l) |(P' JV «> S )= T R (r) y {2) = h { 2 q) /2 (30) 

where p = 1, N p . 

The governing equations at the interior and boundary cells form a system of 48 N p N g linear 
algebraic equations in the unknowns Their solution determines the temperature distribution 

within the FG composite that is subjected to the specified boundary conditions. The final form of 
this system of equations is symbolically expressed as 


kT = t (31) 

where the structural thermal conductivity matrix k contains information on the geometry and ther- 
mal conductivities of the individual subcells (a/37) in the N p N q cells spanning the r and 6 directions; 
the thermal coefficient vector T contains the unknown coefficients that describe the thermal field in 
each subcell, i.e., T = [T^ 11 ^, where = [Tjooo) > 100 ) > TjoiO) > 7j 2 oo) , T(o 2 o) , T(qo 2 ) ]pg^ 

and the thermal force vector t contains information on the thermal boundary conditions. 


NASA/CR— 2001-210967 


8 



4. MECHANICAL ANALYSIS 


4.1. Basic Mechanical Equations. The mechanical dynamic equations of motion in cylin- 
drical coordinates ( r,8,z ) must be fulfilled within each (a/Jy) subcell of the (p, q, s)th cell. These 
equations are given by 


di 7 , 


(<*0y) 


cM a ) 


+ 


da 


(a/3-y) 

Or 


dyW 


+ 




+ 


dzW R( Q 0 y) + f («) 


[(7 


(of07) _ _ p (aP 7) 


'00 


«9 2 u^ 7) 

dt 2 


(32) 


<9r( Q ) dy(P) d zfa} 

dajf^ dat^ dat 0l) 

<9f( a ) <)p(0) <9zC) 


2 

R{ a Pl) -f- f ( Q ) 

1 

R( a Pi) -|- f(a) 


(a/?7) _ „(«07) 
ff rD — P 

n (aPy) _ o (<*07) 

17 rz r 


dt 2 

dht^ 

dt 2 


(33) 

(34) 


where u^ 7 \ u^ 7 \ ui Q/?7) are the subcell displacement components, (ij = r,6,z) are the 

stress components, p- 01 ^ is the mass density and t is time. 

The components of the stress tensor, assuming that the material occupying the subcell (o/3y) 
of the (p, q, s)th cell is either elastic orthotropic or inelastic isotropic, are given by 


a {ocPi) _ c (a07> r e («07) _ e Hafilh _ p(“07> j'(ap-f) ( 35 ) 

ij ijfcl L fcl Ac/ J 

where = r,0,z , c^ 7) are the elements of the elastic stiffness tensor, e^ 7) and e' (Q/?7) are 

the total strain and the inelastic strain components, is the temperature, and T^ 7 ' are the 

elements of the thermal stress tensor which is the product of stiffness and the thermal expansion 
coefficients tensors. In this report, we consider either elastic orthotropic materials or inelastic 
materials which are isotropic in both elastic and inelastic domains. Hence, the above constitutive 
relations (35) reduce to 


= c^ 7) 4f 7) - 2p (a/?7) e-j"^ 7) - <r£ (Q/J7) (36) 

where p(“^ 7 ) is the elastic shear modulus of the material filling the given subcell (0/^7) , and the 
term stands for the thermal contribution . 

4.2. Traction continuity conditions. The continuity of tractions between adjacent subcells 
within the generic cell (p, q, s ) is fulfilled by requiring 


-.007) i(P>7> s ) _ -.007) I (P.9.s) 

ri 'f(l) = d( p >/2 ri 'f( 2 )=-4 P) /2 


(37) 


(0-17) .(p.g.s) _ (027) i(p,9,s) 

'y(l)=/i ( 1 ,) /2 0i y< 2 )=-4 ,) / 2 

(q 01) i(p,g,s) _ (a02) i(p,g,s) 

°zt l 5 <l)=i 1 /2 « lz( 2 )=-i 2 /2 

where i — r,9, z. 


(38) 

(39) 
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In addition to the above continuity conditions within the generic cell, the traction continuity 
at the interfaces between neighboring cells are fulfilled by satisfying 


Mlh) |(p+i,M 


_ cr (2/3'y) |(p,g,5) 


'fW=-<^ ,+1) /2 Ti 'f(2)=4 p) /2 


(ol 7 ) i(p,g+l,s) 


Qi 


a 


(i)=-4« +1 >/2 ° ei l yW=h[ v) /2 

(OC01) | (p,7,s+ 1) _ (a/32) i(p,<j,s) 

(!)=— /i/2 zi 'z^—l 2/2 


(a2j) Ap,q,s) 


(40) 

(41) 

(42) 


4.3. Displacement Continuity Conditions. Similar to the traction continuity conditions 
described above, the following displacement continuity conditions must be satisfied at the interfaces 
within a generic cell (p,q,s) and its neighboring cells. 


„(107) |(M.«) _ ,.(2/37) |(P.?.s) 

'f<i)= 4 p >/ 2 “ U l r( 2 )=- 4'»/2 


s) 


(al 7 ) |(p,9.«) _ ..( 027 ) |(p.9.s; 

u V‘)=^)/ 2 - u '^)= 

(a/31) )(P.9. S ) _ u (a/32) |(p.9.«) 


-4 ?) /2 


u 


l*(i>*i 1 /2~ u U(2) = -l 2 /2 
where denotes the displacement vector in subcell (a/?7), and 


u 


(iPy) i(p+ !,<?,«) __ u (2/?7) i(p><?> 5 ) 




f(2)=4 p) /2 


U 


(al 7 ) |(Pi9+l> s ) 


'y(l)=:-^(« +1 >/2 


_ u (a2 7 ) |(p,9,s) 


|U 

'y(2)=4»>/2 


(a/31) |(P,9,«+1) _ (a/32) |(p.9.«) 

U l*Ci)— Jj/2” U lz(2)=( 2 /2 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


4.4. Boundary Conditions. The final set of conditions that the solution for the displacement 
field must satisfy are the boundary conditions at the top and bottom, and left and right surfaces. 
For example, the tractions in cells (l,g, s) at the bottom surface r = Rq must be equal to the 
applied time-dependent surface loads, 




(49) 


(50) 


where q = and fBr(8,t),fBe(8,t) specify the form of these time-dependent loading func- 

tions. At the top surface r = R\ 


^ Kq».» /a = mm 


(51) 

(52) 
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where q = 1, N q . 

Similarly, the tractions in cells (p, 1 y s) at the left surface 0 = 0 must be equal to the applied 
time-dependent surface loads, 


a 


(aly) 

Or 


|(p.l.s) _ 

'yW^-h.^12 


fbr{r,t) 


(53) 


a 


(“It) i(p,M) _ 
99 1 5(i)=_h ( 1 1) /2 


ho(r,t) 


(54) 


where p = Similar boundary conditions hold at the right surface 9 = 0. In the case of 

prescribed time-dependent displacements (say), they are: 


u: 


(ckly) | (P)Nq t s 


5(2) s 




..(“It) | (p,N q ,s) _ r , .N 


ym=h'f q 72 

where the time-dependent loading functions are denoted by fR r {r,t) and fne(r,t). 


(55) 

(56) 


4.5. Mechanical Field Expansion. The time-dependent displacement components is repre- 
sented in each subcell by a quadratic expansion in the local coordinates f( Q \ y^\ z™ as follows: 


u 


(c* 0 t) 


- w {aPl) 4- f<«)w( Q/5 ^ 4 . ,;(4)t,t/(«^t) 

- ”i(ooo) + r ^ 1(100) + V ^1(010) 

i j(p) 2 i 

+ |(3f<»> 2 - + i(3S OT2 



„(“07) _ 

“0 — 


+ 


”2(000) + ' ’'’winm > y 


2 ( 100 ) 


2 ( 010 ) 


- r(p)2 i mv*/- 1 /* 

i(3f<*> 2 - - -$-)< S! + i(3f™ - i)W. 


ti 


(9)2 


ll. 


4 / KV 2(200) ^ 2 


2 ( 020 ) 


(a/37) 

' 2 ( 002 ) 


(58) 


4 aM = ^tS) + ( s9 > 

where the unknown coefficients (i = 1,2,3), which depend on time, are determined 

from the fulfillment of the governing equations, the interfacial traction and displacement continuity 
conditions, and the applied loading conditions. Note that there are 112 unknowns in eqns (57)-(59) 
which necessitates the establishment of 112 relations for the determination of these unknowns. 

It should be noted that the z- component of the displacement field, eqn (59), does not contain 
linear terms in the local coordinates and This follows from the assumed periodicity in 
the axial direction and symmetry with respect to z ^ = 0 (7 = 1,2). Further, the presence of the 
constant term W^ 0 qq| in eqn (59), which represents subcell center axial displacement, produces 
uniform composite strain e zz upon application of a partial homogenization scheme in the periodic 
direction described in Ref. [5]. This partial homogenization, which couples the present higher order 
theory and an RVE-based theory, leads to an overall behavior of a composite, functionally graded 
in the r and 9 directions, which can be described as a generalized plane strain in the periodic axial 
direction. 
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In the perfectly elastic case, the above quadratic displacement expansions (57)-(59) produce 
linear variations in strains and stresses at each point within a given subcell. In the presence of 
inelastic effects, however, a linear strain field generated by the above expansion does not imply 
the linearity of the stress field due to the path-dependent deformation. Thus the displacement 
field microvariables must depend implicitly on the plastic strain distributions, giving rise to a 
higher-order stress field than the linear strain field generated from the assumed displacement field 
representation. In the presence of inelastic effects, this higher-order stress field is represented by 
a higher order Legendre polynomial expansion in the local coordinates. Therefore, the strain field 
generated from the assumed displacement field, and the resulting mechanical and thermal stress 
fields, must be expressed in terms of Legendre polynomials as 

CO oo oo 

4f rt = EEE )/>„««”) («) 

1=0 m =0 n =0 
oo oo oo 

<r\f y) = EEE (61) 

J =0 m =0 n =0 
oo oo oo 

= w&vu&mw (62) 

{=0 m ~ 0 n = 0 

where i,j = r,6,z, A imn = x/(2 1 + l)(2m + l)(2n + 1), and the nondimensionalized variables £,’s 
defined in the interval — 1 < < 1, are given in terms of the local subcell coordinates as Cf a ^ = 

f(“)/(d Q /2), Cf = ^/(^/2), Ci 7) = * T) /(lr/2). 

For the given displacement field representation the upper limits on the summation in the strain 
expansion (60) becomes 1; while for a quadratic temperature distribution, the upper limits in the 
thermal stress expansion (62) become 2. Alternatively, the upper limits on the summations in the 
stress expansion (61) are chosen so that an accurate representation of the stress fields is obtained 
within each subcell, which depends on the amount of plastic flow. The coefficients m n ^ T ij(i l m,n) 
the a bo ve expansions are determined as described below. 

The strain coefficients e^ Q ,f 7) . are explicitly determined in terms of the displacement field mi- 

1 7 ( i .J71.71 ) 

(a/3y) (&0y) 

crovariables using orthogonal properties of Legendre polynomials. For example e rr ^ 0 0 0 ^ = 

The complete set of nonzero strain coefficients is given in the Appendix. 

Similarly, the thermal stress coefficients can be expressed in terms of the temperature 

field microvariables Tfj^y For example T'^oAO) = The com P^ cte set °f nonzero 

thermal stress coefficients is also given in the Appendix. 

The stress coefficients are expressed in terms of strain coefficients, the thermal stress 

coefficients, and the unknown inelastic strain distributions by first substituting the Legendre poly- 
nomial representations for e ^ 0l) , and ^(a/h) into the constitutive equations (36) and then 

utilizing the orthogonality of Legendre polynomials. This yields 


(a 0 ~t) _AaPi)J°‘Pi) _T{aPl) _ 

T -ijko e /co(i ,771,71} 1 ij(l,m y n) ij (l ,171,71) 


(63) 


The Awf 7 ' 1 x terms depend on the inelastic strain distributions calculated in the following 


manner: 
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p (a/? 7 ) _ .. A-l mn 

K y(i,m,n) — ^((a/ 37 ) ^ 



4- Q/37) P; (cl°° ) p m (cf } (d 7) ) d cl°° dc? } d ci 7) 


(64) 


4.6. Stress Moments. Let us define the following stress quantities in subcell (ct/?7) of cell 

(p,q,s) 




Jp.g.s) 


f f 0 f (r^ y (t/^) m (2^ 7 ^) n (r|° ! ^ 7 ^dr^ a ^ dy^ dz^ 

-d^ ) /2J-h^ ) /2J-ly/2 ] 


(65) 


where i,j = r, z. Note that, in particular, the zero-order quantities S^ Q 0 0) represent the average 
stresses in the subcell. 

Explicit evaluation of eqn (65) yields the following expressions for the normal stress quantities 
(superscripts (p, q, s ) have been omitted) 


0 ( 0 / 37 ) 

‘-’rr( 0 , 0 , 0 ) 




-11 


1 ( 100 ) 


+ c! 


12 


] [W, 


2 ( 010 ) + £>(a/ 3 7 ) t ' K 1 ( 000 ) 


W\ 


(«/? 7) 


] + C 


(a/37) 
13 


Ty(«/37) 
^ 3 ( 001 ) 


'p(ap'y)rp{&P'y) __ jA&Pl) 

1 rr ^ ( 000 ) ‘‘Vr ( 0 , 0 , 0 ) 


( 66 ) 


with similar expressions for 


°00(O,O,O) 


and 

ana ^(0,0,0) ■ 


0 ( 0 / 37 ) 

*^rr(l, 0 , 0 ) 


dl 


,(o^7)tt/(o/37) 


-11 


1 ( 200 ) + 


J.O‘0l)p<r( a 0l) 

12RW-1) 12 1(100) 



da n (a/37) 

2^3 rr(1 ’°’ 0) 


(67) 


with similar expressions for 


q( a @y) 

° 69 ( 1 , 0 , 0 ) 


and 9 (q/37) 
ana o zz(l, 0 , 0 )' 


q(a/ 3 y) 

°rr( 0 ,l, 0 ) 


fys (a/ 3 7 ) 

T 12 


w(«^) 

vv 2(020) 


+ 


r (a0y) w (<x01l) 

l2R( a Py) ° 12 v 1 ( 010 ) 


(*0i)rAaPi) _ JV_d(q/9 7 ) 
12 irr ( 010 ) 2 v ^ /l rr( 0 ,l, 0 ) 


( 68 ) 


with similar expressions for and 


o( Q / 37 ) _ ^a r o(a/ 37 ) 

a rr( 2 , 0 , 0 ) — ^ ^"(O.O.O) 


o(a/37) _ r o(a/37) 

^rr (0,2,0) ^ *■ rr(0,0,0) 

q( a Py) _ ix 10 ( 0 / 37 ) 
J rr(0,0,2) ~ j 2 [ ^(0,0,0) 


2 (a/ 37 ) 9 

a ( c 12 71/(0/37) -p(a/ 37 ) rpfaP-y ) \ £_ rf(apy) 

0 l 7R(a/37) ”1(200) 1 rr J (200) / ^"^(2,0 


, 0)1 


/ cjf 7) TT/ (a/ 3 7 ) _ r ( a py)r r (ap-,) ) 2 _ n(afiy) 

0 ^R(aPj) ^ 1(020) 1 (020) ' ^^(0,2 


(0,2, 0)1 


(a/ 37 ) 9 

t C 12 tt/(o/ 37) _ r (a/37) T (a^7)\ f_ rp(a/37) 

V 1 1 rr ^^0021 / S-rfO-O. 


with similar expressions for the other components. 

Similarly, the explicit expressions for the shear stress quantities are 


(69) 

(70) 

(71) 
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c(“/3 7) 
c ’r0{O,O,O) 



1 

R(a0 7) 



jy( a 0y) 

■ fl r0(O,O,O) 


q(a0y) 

°r0(l,O,O) 


c(a/J 7 ) 

°r0(O,l,O) 


^o „( a 0y) \l{ri a Pl) I W (o/3 7 ), _ p (a/3 7) 

4 44 1 2(200) 3^(o/3 7 ) ^2(100)1 2V3 r ° (1 ’°’ 0) 


4 ° 44 


(«^7) ^(“07) 


1 


1(020) 3^ 


w( Q ^7)l 

:/3 7 ) ^ 2(010)1 


fyj n{a0l) 

2 ^3 r0(o - 1,o) 


q(a0y) _ djrq(alh) ffa Jp0i) Uf( a 0l) qi&Pl) 1 

r0(2,O,O) _ J 2 l °r0(O,O,O) 10/?(^7) 44 ^2(200) ^ n r0(2,0,0 ) J 

q{a(3y) _ Z'qXPy) (a/3 7 ) ti/(o/3 7 ) 2_„(a0y) 1 

°r5(0,2,0) “ 12 L°r0(O,O,O) 10i?(°/97) 44 2(020) / V0(O,2,O)J 

q{a0y) _ ix fo(o/3 7 ) {7 >Wm(aW ft{«0 7) , 

°r0(O,O,2) “ 12 l°r0(O,O,O) 10 fl(a/3 7 ) C 44 ^2(002) ^^(0,0,2)! 

q{a0y) _ ^7 „(o/37)rj/(a/37) _ _Jt_ d (<*07) 

°rz(0,0,l) ~ 4^44 KK 1(002) 2 ^ rz ( 0 >°> 1 ) 

q(o/3 7 ) _ ^7 (a/3 7 ),,,-(a/3 7 ) _ ^7 Tf{a0y) 

°02(O,O,1) — 4 c 44 ^ 2(002) 2 v / 3' n ' 0z ( O ' O > 1 ) 


(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 


4.7. Zero-Moments of the Equations of Motion. By integrating the three equations of 
motion (32)-(34) over the subcell (a/?7) we obtain, in conjunction with the above displacement 
expansions (57)- (59), 


T (a0y) 


+ J, 


(a/37) 


+ K 


(a0y) 


rr(0,6,0) ‘ r *"0r(O,O,O) ^ zr (0,0,0) + ^j.(a/3 7 ) 


7(0^7) , f 

r0(O,O,O) + -"00(0,0,0) 


,{<x0y) jA<*@y) 1 

1 /V n n\ i~ 


fc( a W o(o/? 7 ) 1 _ ^(a/37) T j/(a/37) 

n D rr(0 1 0,0) ^(0A0)J“^ ^1(000) 

,c( a W — 0 ( Q /57)w( a ^) 
^rflro.n.m ~ t * yv 2(000) 


z0(0fifi) R {ai3 7 ) °r0(O,O,O) 


r(«^7) 


+ Ji 


-(a/3 7 ) 


~(a/? 7 ) , 1 

1 nm T 


7(0/37) 


"72(0,0,0) ' r -"02(0,0,0) T ^22(0,0,0) ^(o/3 7 ) O r2(0,0,0) 

where the following interfacial traction integrals have been defined: 


-0 


®0) = 2^-4)” r' 2 f"' 2 rifAy) + (-i)" + v<f i) (-|)i<is OT ^w 

; ^(a/? 7 ) 2 J-hp/2J-ly/2 J 2 ^ 2 

«.o, ~ 7 >- JZ Ov*® + 

- zA r fzizy^ + 


(80) 

(81) 

(82) 

(83) 

(84) 

(85) 


where y = r, 0, z\ n = 0 or 1, and a^ 7 \±d a /2), a^^ y \±hp/2), <r*“^(±Z 7 / 2), stand for the 
interfac.ial stresses at ±d a /2, ±hp/2, ±1 7 / 2, respectively. 
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4.8. First-Moments of the Equations of Motion. By multiplying the three equations of 
motion (32)-(34) by f^ a \ and integrating over the subcell volume by parts, we obtain: 


Aapy) q(aPy) . 1 r 0(0/37) q(a 0 l) ) _ Ja 0 y) ^o pr/(«/37) 

i rr(l, 0 , 0 ) °rr( 0 , 0 , 0 ) ^ R{apy) l°rr(l, 0 , 0 ) ° 00 ( 1 , 0 , 0)1 ~ P 42 1 ( 1 °°) 

( 86 ) 

T (aPl) qi.aPl) , 2 q(apy) _ Ja 0 y) ^q rx/i^Pl) 

r0( 1,0,0) °r0(O,O,O) * r R(apy) °r0(l,O,O) ~ P 42 2(100) 

(87) 

A a Pl) oi a 0 l) 1 1 o(q/ 37 ) _ (1 

i r 2 (l, 0 , 0 ) °r 2 ( 0 , 0 , 0 ) ‘ r ^>(0/37) *-Vz(l, 0 , 0 ) u 

( 88 ) 

By multiplying the three equations of motion (32)-(34) by t/^, and integrating over the subcell 
volume by parts, we obtain: 

1 h 2 

r(a/ 3 y) q(aP 7) . 1 o( a ^7) 1 _ 

J 0 r(O,l,O) _ J r0(O,O,O) ^ ^(q/3 7 ) l°rr(0,l,0) °00(O,1,O)J ~ P 42 1(010) 

(89) 

j( a Pl) q(aPy) . 2 0 ( 0 / 37 ) _ -(a/37) ^0 \jA a Pl) 

,y 00(O,l,O) ^(/(O.O.O) ^ #(a/3 7 ) °r0(O,l,O) “ P 12 *2(010) 

(90) 

r(o/37) r.(o/37) , 1 o(“/37) _ n 

J ^(0 t l,0) ° 0 *(O,O,O) + ^(a/3 7 ) °^(0,1,0) “ U 

(91) 

By multiplying the three equations of motion (32)-(34) by z^\ and integrating over 
volume by parts, we obtain: 

the subcell 


K 


(ocf 3 y) 

r(0,0,l) 



1 r g{aPl) 
R(aP 7 ) l a rr( 0 , 0 ,l) 


q( a P l) 
°00(O,O,1) 


] = 0 


jr(a(3y) _ q(apy) £ q(apl) 

n z0( 0,0,1) D 9z( 0,0,0) ' 1 ~ R(c*l3y) °r6(0,0,l) 


= 0 


Tf(aPy) q{ap 7) 1 qi^Pl) 

^ 22 ( 0 , 0 , 1 ) *" > 22 ( 0 , 0 , 0 ) i “ R(aPy) °r 2 ( 0 , 0 ,l) 


P j 2 ’^ 3 ( 001 ) 


(92) 

(93) 

(94) 


4.9. Second-Moments of the Equations of Motion. By multiplying the three equations of 
motion (32)-(34) by f^ 2 , and integrating over the subcell volume by parts, we obtain: 


^o Aa0l) 

4 - 1 rr( 0 , 0 , 0 ) 


+ 


+ 


t(o/3 7 ) 

12 t J 0r(O,O,O) 

1 

R{<*0 7 ) 


+ K 


( 0/3 7 ) 
2 r( 0 , 0 , 0 ) 


o o( Q h7) 
^°rr(l,0,0) 


[ c( Q Py) _ C'( q /37) 1 _ JocPy )^[ph(°'/ 3 7) , £qpb( a ^7)l 

1*^(2, 0,0) 1 ' 5 ( 9 e( 2 , 0 , 0 )i ^ j^'-^qOOO) w ^i(200)i 


dl 


(95) 


4 Mpy) 

4 J r0(O,O,O) " t ‘ 


+ 


€ 

12 ^ee( 0,0,0) 

2 q(a/3y) 


! f A 0 Pi) 
il J 06 

-S K _ 


+ K 


(apy) 


20(0,0,0) 

R(a0l) *“V0(2,O,O) - ^ 


1 _ O <?( Q ^7) 

1 Z °r0(l,O,O) 


(q/? 7 ) ((a itV/(o/ 37) , (^opr/O^h 

42 ^ 2 ( 000 ) + 49 2(200)J 


^q T (a0y) , ^q r t(o/37 ) 

4 i r 2 (o,o,o) 42 ^^z(o.o.o) 


+ K 


(aPy) 1 

22(0,0,0)1 


-25. 


(q Pi) 


r 2 (l, 0 , 0 ) + p£(aPy) Lj rz(2,0,0) 


Aapy) 


= 0 


(96) 

(97) 
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By multiplying the three equations of motion (32)-(34) by and integrating over the subcell 
volume by parts, we obtain: 


j(<*Pl) , ^/3 [ , ffiaPl) 

4 J 0 r(o,o,o) ' r 12 [ rr (°’ 0 ’°) ^ ' 


zr(0,0,0)J 


9 q(^Pl) 
Z °r0(O,l,O) 


+ 


1 




(a/?7) 


J^(a0~f) l^rr(0,2,0) ^5(0, 2,0) 


q( a Pi) | _ (a/3 7 )^f 14 /( 0 / 37 ) . ^hpp(“/ 5 'r)i 

C, W0.2.0)I — P 12 1 1(000) ^ IQ ' ^ 1(020) J 


(98) 


^ 0 j(<*P 7) , r r(«/37) , K {aP~l) 1 _ 9 q{a0l) 

4 J e 0 (o,o,o) + 12 ^*^( 0 , 0 , 0 ) -“z 0 (o,o,o)J z ° 00 (o,i,o) 

4 . 2 c(“/^7) _ (a/3 7 )^hrrir(a/?7) , ^pj/(o/?7)i 

+ £(a/3 7 ) O r0(O,2,O) “ P 1 2 l ^ 2 ( 00 °) + 10 2 (° 2 °) J 


^0 j(aPif) . hfi I AaPi) . jsiapi) 1 _ 9 q{ a (h) , ^ q(o/37) _ 

4 * / 0z(O,O,O) + 12 l i rz(0,0,0) -+ ' -^zz(0,0,0)J ZO 0z (0,1,0) + ^(a/3 7 ) °rz(0,2,0) 

» i-t, ^ /QO\ f O. A\ Iatt nn/4 1 rtrr»> ri f in rr /7' 


(99) 

( 100 ) 


** ' ■ ■ ' j. ^ ' j . t ' ■ 

By multiplying the three equations of motion (32)-(34) by zW 2 , and integrating over the subcell 
volume by parts, we obtain: 


^7 T^iaP'r) , 4 r Aafry) , Aa0i) 1 _ 9 o(nPl) 

J ft Z7-(0,0,0) + i 2 t i rr( 0 , 0 , 0 ) ■** J 6»r(0,0,0)J ZJ rz(0,0,l) 

1 1 r q( a Pl) q( a Pl) 1 _ Japi 

* TZTZoZX [J rr rn n o\ ‘-’flam n 91 J ~ H 


1 / 2 7 2 

1 r q( a Pl) q{ a P l) ] _ Xa/ 37 ) J \\xA a P^) , T_\ir{ a Pl)] 

^(a/hO o rr(0,0,2) °00(O,O,2)J ~ P 12 1 ^000) ^ 10 KK l(OO2)J 


( 101 ) 


& 

4 A z0(o,o,o) 




+ 

+ 


ll 


rH“P7/ ■ t(“/?7) 1 nqK<xpi) 

_ l r0(O,O,O) J ^(0,0,0)) Z pz(0,0,1) 

2 q(<*Pl) _ Ja/37) ix fw( Q ^) + ix ppW 5 '*)} 

^(ft/37) °r0(O,O,2) ~ P 12 1^2(000) + 1Q y * 2(002)) 




/ 2 

S ^( 0 ^ 7 ) , '7 [ 7-(oP7) , 7 l“P7l 1 _ 9 q(,“P7/ , * M“P7/ _ q 

^(0,0,0) ■*" 22 ^ 2 ( 0 , 0 , 0 ) ^ J ez(0,0,0)J ZJ zz(0,0,l) ^ R(<xfa) °rz(0,0,2) 




r(a/?7) 


z(“h7) 


(a/37) 


j(<*/37) _ 


( 102 ) 

(103) 


4.10. Interfacial Traction Integrals. Prom the above 21 relations (80)-(82) and (86)-(103) it 
is possible to establish expressions for the interfacial traction integrals /h( 0 , 0 , 0 )> ^(o, 0 , 0 )’ ^ z;(o,o,o)’ 
7 = r.O.z, as follows. — 

J y ’ ''j t # (ctQ'y) 

Substitution of eqn (80) into eqn (95) yields the following expression for I^/o.O.O) ' 


t(<*P 7) 
i rr( 0 , 0 , 0 ) 


12 MaP'f) , 1 r q{a0 7 ) q{a0l) 1 

^2 *^rr(l,0,0) ' 2J?( Q / 5 7) '° rr ( 0 , 0 , 0 ) 0 00 ( 0 , 0 , 0 )) 


£ f c(<>/57) 

^2^(a/37) )°rr-(2,0,0) 


q ( q / 57 ) 

o 00(2,O,O) 


hi 


1 4 . £,( 0 / 37 ) 1 ^ 01 )^( 0 / 37 ) 
J + P an 1(200) 


20 


( 104 ) 


Similarly, substitution of eqn (81) into eqn (96) yields 
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12 Q (ap-y) 


1 Q (o/3 7 ) 


V0(O,O,O) “ ff 70(1,0,0) _r ^(a/?7) ^ 70(0, 0,0) 

12 cM 7) | „(a/?7) II/(“/ 3 T') 

rf2 ^(a/37) ^0(2, 0,0) + P 20 2 ( 200 > 


Substitution of eqn (82) into eqn (97) yields 


j( a Py) 

•Vz(0,0,0) “ U 


(afiy) . 


Substitution of eqn (80) into eqn (98) yields the following expression for 

T (a0y) _ 12 Aafiy) . 1 [AW _ c(a0y) 1 

J 6r (0,0,0) ~ ^2 *70(0,1,0) 2R( a ^ ^ rr ( 0 . 0 . 0 ) 00(0,0,0)1 

6 r qi a Py) e(o/37) 1 1 n (aW^w( a W 

” 1 "(0,2,0) *00(0,2,0) J ' P 20 1(020) 


Similarly, substitution of eqn (81) into eqn (99) yields 


j(<*0y) _ 12 q( a /37) , 

^00(0,0,0) — u2 ^00(0,1,0) ' r 


1 0(0/37) 

R(afiy) 70(0,0,0) 


12 ^-,( 0 / 37 ) ■ Jafiy) hfi tj/(°/37) 

hiRiafh) ^e(0,2,0) + p 20 2 (° 2 °) 

Substitution of eqn (82) into eqn (100) yields 

t(o/ 37) _ n 

J 02(0,0,0) — u 

f q/^ i "v1 

Substitution of eqn (80) into eqn (101) yields the following expression for K zr{Q Qfi) : 

Ts ( a P 7) _ 12 0(0/37) I 1 [ 0(0/37) _ 0(0/37) | 

A 27(0,0,0) “ 12 * 72 ( 0 , 0,1) + 2/?(o/37) [ °rr(0,0,0) °00(0, 0,0)1 


6 r 0 ( 0 / 37 ) 0 ( 0 / 37 ) 1 _i_ Ja0y) ^7 tt/(o/37) 

(S^ 1*77(0,0,2) *00(0,0,2)1 + P 20 1 (° 02 ) 


Similarly, substitution of eqn (81) into eqn (102) yields 

( o / 3 7 ) _ 12^(0/37) , 1 0(0/37) 

21 20(0,0,0) - ^2 *02(0,0,1) ' J^(a0-y) 70(0,0,0) 

12 0 ( 0 / 37 ) , „(o/3 7 )ix br/(“^7) 

l2R(a/3y) *70(0, 0,2) + P 20 ^ 2 (°° 2 ) 

Substitution of eqn (82) into eqn (103) yields 


jAaPy) _ n 

22(0,0,0) — u 


(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 

(HI) 

( 112 ) 
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It should be noted that the remaining nine expressions for the interfacial traction integrals 

7 rj(i AO) > J dj(o},o) ' K ^(ol,i) 0 = r ' e ' z ) have steady been given by eqns (86)-(94). 

It is convenient to define the following quantities 


A<*Pj) _ A^P 7) Ju0y) 4 iiA°‘P'l) 

Vr(0,0,0) ~ i rr(0,0,0) P 2 0 M 2 00) 

(113) 

A a 0 7) _ rfapy) n (a0y) ^/(apy) 

r8( 0,0,0) ~ r9(0,0,0) P 20 KK2 ( 20 °) 

(114) 

A a P 7) _ Aa0~t) Ja0y) 4 uA a 0'r) 

Vr(l,0,0) ~ i rr(l,0,0) " Kl l( 100 ) 

(115) 

.(a/97) _ r(Q0 7) Japy) 4 yrr(apy) 

r9(l, 0,0) — J r9(l,0,0) P 12 2(100) 

(116) 

Aa0y) _ Aa0y) Ja0y) 4 yy(apy) 

J6r(0,0,0) — ■ y 9r(0,0,0) P 2 0 ^(020) 

(117) 

(a/?7) _ Aa0 7) (a/3 7 ) 4 

J 88(0,0,0) ~ J 66(0,0,0) P 20 2 (° 2 °) 

(118) 

.(a^7) _ j(a0y) Jotp-y) 4 rV^(a/3 7 ) 

“'9r(0,l,0) “9r(0,l,0) P ^^H 010 ) 

(119) 

.(q^ 7) _ Aa0y) Ja0y) 4 At/i a Pl) 

•>06(0,1,0) ~ J 88( 0,1,0) P 12 *^2(010) 

(120) 

,(a0j) _ K (a0l) (a0y) ^7 w( Q ^7) 

K 2r(0,0,0) ~ -^zrfO.O.O) P 20 M 002 ) 

(121) 

j2 

i(a0y) Ts( a 0~t) (a/37) S w( a ^7) 

^29(0,0,0) — 29(0,0,0) P 20 KK2 (° 02 ) 

(122) 

, (a/37) _ js(apy) >/37)iw(“W 

^22(0,0,1) — -^22(0,0,1) P 12 K *3(001) 

(123) 


It should be noted that the above eleven new quantities (113)-(123) are given solely in terms of 
the stresses eqns (66)- (79), and therefore involve no time derivatives of the unknown 

displacement coefficients 

4.11. Volume Averages of the Equations of Motion. Once the interfacial traction integrals 
(83)“ (85) have been established, we can readily express the volume averages of the three equations 
of motion, given by eqns (80)-(82), over the (a/?q) subcell of the (p,g, s)th cell. As it can be readily 
observed, eqn (82) is trivially satisfied, while eqns (80)- (81) are, respectively, given as follows: 




1 ( 000 ) 


4 

20 


(a/?7) _ 4 r irjaP-l) 
1(200) Of) 1(020) 


l 2 

_L 

20 ' 


_ IwWi) 

"" '* 1(002) 


A a 0i) 

l TT 


ap-r) , A^py) , T.(<*P1) , r c( a Pl) _ q\<xpi} 

r(0,0,0) + .?9r(0,0,0) ' zr( 0,0,0) + R(a0y) Prr^.O.O) ,:> 99(0,0,0) 


A a Pi) 


(apy) 


gi'xfi'r) 


(124) 
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Ja0~i) |ir/( Q 07) I 

p i M/ 2(O0O) 20 ^ / 2( 200 ) _ 20 ^C 020 ) 20 VV2 (° O2)J 

_ A<*Pi) ,A a Py) i i.( a 0~<) | ^ g,(a/?7) ('125') 

— Ve(o,o,o) + ^00(o,o,o) + K z0(o,o,o) + j^ a p T ) »-0(o,o,o) v ; 

These two equations form 16 out of the 112 relations needed for the determination of the second 
time derivatives of the unknown displacement coefficients W^ lmn y 

4.12. Imposition of the Traction Continuity Conditions. The traction continuity condi- 
tions are imposed on an average basis at the subcell and cell interfaces. These conditions imply 
existence of certain relationships between the aforementioned interfacial traction integrals as de- 
scribed next. 

Interfacial continuity of tractions in the radial direction. Let us define the following 
two new quantities: 


j(P) ; 


(126) 
(127) 

rj ■ >j 'fW=d;T7 2 ' J /z 

with j = r , 0 , z. 

Substituting the interfacial traction continuity conditions (37) and (40) between cells and sub- 
cells in the r-direction, we obtain, respectively, 


.(olP 7 ) 

1 (p,g. 3 ) 


| (p,9,«) 

rj 

'f(“)=4 p) /2 

°rj 


.(a£7) 

|(Pi9i s ) 

+a£f 7) 

rj 


rj 

l f=(“)=d<f ) /2 

f ( a )=— 1 


F APy) I (p, <?,«)_ J2M) I (?,<?>«) -o- ( 2/37) |( p_1 > 9 >/) 

r Tj I rj 'f{2)=-4 p) /2 rj l f(2)=4 p - 1) /2 

r (iPl) j (p,q,s) J2P-t) |(p,?, S ) |(P- 1 '9.s) 

G rj I ~ a rj l f (2) = _ d (p)/ 2 +a rj /I 


(128) 

(129) 


By addition and subtraction of equal quantities to and from the last two equations it can be easily 
verified that 

2F (1/?7) |(m-*)= + Gj.f 7) ] (M ' s) - {Flf l] + 7) ] (p ~ 1,<? ' s) (130) 

2 qAW [~F^ y) + G ( ^ 3l) } (p ' q ' s) + [Fj:f l) + G^ 7) ]( p_, ’ 9,s) (131) 

Then employing the definition for e 0 n ( 83 )» we obtain the corresponding relations: 

r( 1^7) |(p,<j,s) _ _y j( 2 P~f) |(p,<}, 3 ) , |(P -1 '9. S ) 

2 rj(0,0,0) I — Zl ''rj(0,0,0) I ^^(o.o.O) ■ 

4 - 7,r (2/?7) |0 vj> 3 ) -Z,r^ 7) Kp- 1 ^'®) (132) 


7) |(p,<7, 3 ) _ Kp.<7. 3 ) [(P" 1 -?, 3 ) 

j(l,0,0) I " 5j rj( 1,0,0) I ' l 


rj(l 


y t{ j(p,-J,s) i 7 r(2/57) |(p-l, 9, 3 

“ ^^(0,0,0) i + Z 8i ri(0 ,0,0) I 

where j = r, 0, since for j = z all quantities vanish. 


(133) 
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The above equations (132)-(133) are applicable for internal cells where p = 2,3, .„,N P . For the 
1st cell p — 1 the continuity of tractions at the interfaces between the subcells of this cell can be 
easily shown to yield the single relation: 

[tf£f 7) + F^)}^) = - F^ ] f' 9 ' s) (134) 

This single relation replaces the above corresponding two relations (132)-(133) at interior cells. 
This relation provides 


r jQfh) i JL t^Pj) 1(1,9, s) _ 7 f_z, 

‘ rj( 0 , 0 , 0 ) + dl i r;(l, 0 , 0)1 “ “1 * ll rj ( 0 , 0 , 0 ) + Z ' 3 - l rj(l 1 0 1 0 )J 


(135) 


where j = r, 9 , since for j = z this equation is identically satisfied. Again, this relation replaces the 
above two equations that are applicable at interior cells. 

Substituting the established expressions for the interfacial traction integrals into relations (132)- 
(133), we obtain the following equations for internal cells when j — r 








20 kk 1(200)J -VIV* 20 1(200}J 

+ UP [20l) %Mt& V ' q ' S) ' Z 


12 1 ( 100 ) J 


= l 


;(1 fH) 

rr(0,0,0) 


|(p, 9 , s ) J_ V A 2 / 3 ^) |(p, 9 ,») 4-Zoi^ 2 ^ |(p-l, 9 ,s) 

I ' ^l*rr( 0 , 0 , 0 ) I +Z/ 2 Vr fO.O .01 I 


rr(0,0,0) 






12 1(100)J 


20 1 ( 20 °) J 


_ AWi) 


= i 

rr(l 

+ Z*ji 


f) \(p><h s ) _ Z l(p^ 5 ) i(p-bq^) 

,0,0) I ^5 Vr(l,0 t 0) I ^ 6 Vr(l,0 3 0) I 


;( 2 0y) i(p,g,s) i(p~ 1,9,5) 

I ** l rr(QA.O) I 


*rr( 0 , 0 , 0 ) I " ' ^ 8 * rr (o, 0 , 0 ) 

For j — 6 , the two relations (132)-(133) become: 


I™] 1 '''’'*’ - 


+ 




+ 


_ ,-(1/37) 

r<J(l 


= t. 


f) |(p,9,s) -Zri^ 2 ^ ifp- 1 .?.® 

,0,0) I "5* r 0( 1,0,0) I z>6 V0(1,O,O) I 


+ Zji 


(2/37) 

r 0 (O,O,O) 


|(P, 9 ,«) _ Z 


9(0, 0,0) 


|(P-! 


-1,9,*) 


(136) 


+ Zslr> m § H'fjP <ww) + z 6 [^4» , i < £)I < p " i,,,,) 




(137) 


5) 


«au l (w) +z ‘d,0) I''’ -1 "' 

Z^l.0) +^ i S(W,0) 


(138) 


(139) 
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Consequently, these four equations (136)-(139) provide additional 16 relations for the second time 
derivatives of the unknown displacement coefficients W^ lmr ^y 

The two relations (132)-(133) are applicable at internal cells p = 2, N p . For the boundary 
cell p = 1 we use the corresponding relations (135). For j = r we obtain 


1 

+ 

2Z 1P ^§W%™ 



— 




— 

22.^0, o, - 1*^,0) + 223^0, 

(140) 

while for j = 0: 




iP 20 ^ 2 ( 200 ) 

+ 

22 1P <«|^W + wSm 



— 




- 

07 AW-y) 2 AWi) , n 7 J 2 ^ l(i>g,s) 

^ Zl V0(O,O,O) ^ Vl9(l,0,0) + Z//3 V0(1, 0,0)1 

(141) 


These are 8 relations that are applicable for the boundary cell p = 1. For this cell however there 
are another 8 relations that express the applied boundary conditions at the surface r = Rq. These 
latter relations will be presented in Section 4.14. Consequently, we have altogether 16 relations for 
the second time derivatives of the unknown displacement coefficients in this boundary cell, 

p = 1, (just like any inner cell p ^ 1). 

Interfacial continuity of tractions in the angular direction. A similar analysis for the 
traction continuity conditions in the ^-direction yields for internal cells q = 2,3,..., N q 


Aaly) |(p,g,s) 

J 9j(0,0,0) I 


— v, r(“ 2 'd i(p.« 7. s ) 


+ VsJE 


'<b( 0 ,o,o) 

(a2-y) 


-YiJZ 


(a2y) |(p, 9 -l,s) 


07 ( 0 , 0 , 0 ) 


0j(O,l,O) 


l(p 


P,g ,s) _Y 4 j( a ?y) |(p>9 — M) 


( 0 , 1 , 0 ) 


|(P 


(142) 


j(al~r) |(p,q,s) _ y t( q2 7) |(p,<J,s) iV |(p,q-l,s) 

J ^(0,l,0) I “ y 5«V, (0 ,l,0) I +^^(0,1,0) I 

— y, T^ 2 ^ |(p.9. s ) 4 -Ya J ( ' a2 ' r ' 1 x |(p.'7-i, 5 ) 

YlJ 9j{ 0,0,0) I +y 8U 0j(ooo) | 


(143) 


where j = r, 0, since for j = z all quantities vanish. The above equations (142)-(143) are applicable 
for internal cells where q = 2, 3, ..., N q . 

For the 1st cell q = 1, the continuity of tractions at the interfaces between the subcells of this 
cell yields the single relation: 


r x(ofl'y) i JL i(p,M) _ o[_y t(" 2 7) _l y, /( a2 'i') ] (P>1» 

[^(0,0,0) + ^0j(O,l,O)J — (0,0,0) + r 3 J ^'(0,l,0)J 

where j — r, 6 , since for j — z this equation is trivially satisfied. 


i«) 


(144) 
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Substituting in (142)- 143) the established expressions for the interfacial traction integrals we 
readily obtain for internal cells the following expressions when j = r: 


+ n[p (tt27) ^g5 f ft p ’ q ’ s) - n[^ 27) l<S i ] (M - 1 ' s) 


12 i(oio) J 


— ,-(“ l T r ) |(p,9,s) , y, „-(«27) |(p,<j,s) , y n1 '(« 2 '») |(p,<J-M) 

•'5r(0,0,0) I + - ri -'5r(0,0,0) I ' 2 ^0r(O,O,O) I 

— Yii^ 012 ^ |(P-9>s) i y.^ 02 '!') j(P,9-l.s) 

1 3-'0r - (O,l,O) I ^ _r4 -'5r (0,1,0) I 

_ |(p, 9 ,s) _y c1 ( a2 7) |(p,9,s) _y, 1 ( Q2 ')') |(p,g-l,s) 

— •'5r(0,l,0) I ' ft- (0,1,0) I I GJ0r(O,l,0) I 

+ - y < 2 4.o) 

For j = 5 the two relations (142)-(143) yield, 


— f 0 ( al 7)^w( Ql ^)l(p»<?,5) _ Y f 0 ( a 27)^2 tt/(Q£27)i(p^ > s ) Yo[^ a27 ^ — 

^ 20 2 (° 2 °) on K *2f020P ^21^ on KK 2f020W 


20 2(020) J 


20 2(020)J 


+ >3[p (a27) ^^2(010)]( p ’ 9,s) - 

_ ? -( al 7) |W ±7,/“^) i(pi 9> s ) xy,^°^ |0,<?-i,s) 

-755(0,0,0) I 1-755(0,0,0) I + J 2-755(0,0,0) I 

_ y ,(o27) |(p.<m) I y l1 ( Qf2 T') |(p,9-l,s) 

^ 3 -^ee>(o,i,o) I ^ Z4 J00(o,i,o) I 




12 2(010) J 


12 2(010) J 


— i( q 1 7) |(p,9,s) _y 1 ( a2 7) |(p,9,s) _y „'( Q2 7) i(p,y-l,s) 

- •'55(0,1,0) I J 5Jw(0,l,0) I ^6^(0, i,0) I 

+ l (w> - y «&$.o) l (P ' , " 1, ’ > < 148 ) 

Consequently, the conditions that the tractions are continuous at the interfaces in the 6 direction 
provide 16 additional relations which are valid for internal cells q = 2,3, . . . , N q . 

For the boundary cells q = 1, we obtain from eqn (144) the following expression when j = r: 


r „(al7) ^1 U/( al 7) , ny Ja2y) ^2 Tj/(a 2 7) , „(al 7 ) Txr( al 7) 

IP 20^1(020) + 2Y 'P 2Q W l{020) + P 6^1(010) 

— V^ a2 ^W^\ 1 (?.!,«) _ f_ ? (“)7) 

R KK 1(010)J l J5r(0,0,0) 


oy i( a2 7) f_ -( Ql 7) 4-2 Vq’i^“ 2 ^ lfp- 1 ' 4 ) 

Z1 W57-(0,0,0) 0,1,0) + Zr3 J5r(0,l > 0)l 
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while for j = 0: 


r / ,( al 7)^lw( al T , l 4. 2Yi n( 0 ‘ 2 ' y ' 1 — 4- 0 ( al 'i') — W^ a1 ^ 

IP 20 ^ 2 (° 2 °) + AYlp 20 ^ 2 ( 020 ) + p 6 2 (° 10 ) 

- = 0) 

- 2rd^ i0) -^5;l,„) + 2y3j% ) 1 , () )] <,, ‘ 1 ’* ) (iso) 

Thus for the boundary cells q — 1 we have 8 relations which together with the 8 imposed boundary 
conditions at the surface of this cell 0 = 0 (which will be discussed in Section 4.14) form the required 
16 relations (just like any internal cell q ^ 1). 

Interfacial continuity of tractions in the axial direction. Let us consider the traction 
continuity conditions in the periodic axial z-direction. Here we define the two quantities 

\(p^s)_ {Otp 7) |0 *,*.*) J*Pl) |(M,a) (1 51 \ 

b zj V - a zj U(7)^ 7 /2 ~ a zj l*C 7 )=-I 7 /2 

r*( a Py) I (p^s)_ n (<x 0 l) 1 (P, 9 ,«) . i(p,g,a) f 152 ^ 

I l*(7)=/ 7 /2 iz(7) = -i 7 /2 V } 

with j = r, 2. It can be readily established that due to periodicity of the stresses between repeating 
cells in the axial direction, the continuity of tractions at the interfaces between the subcells of the 
cell yields 

p(ot/ 3 l) j(p,g,s) __ _p( a P 2 ) |(jV7,s) 

Q(afil) |(p t g,5) = q{<*P2) |(p,< 7 ,s) 

The above first equation (153) gives the following nontrivial relations 

0,0,0) + 0,0,0) J u 

for j = r and 0, while for j = z we obtain from eqn (154) 

f Tfiotfil) tv-(q/^ 2) i (p,q,s) _ n 

l ^ zzi 0,0,1) ^**(0,0,1)! “ u 

Consequently, the following relations can be readily established in the cell (p, 9, s ) 

(a/ 3 l) ^1 rzr(at^l) , (a/ 52 ) _^tVt(^/ 32 ) . . (a/ 31 ) . v (a/ 32 ) 

Z 9 20^1(002)+^ 20 V ^ 1 (° 02 ) - Zl ^r(0,0,0) i2 Sr(0,0,l 

for j = r, and 

(a/ 31 ) ^1 w( a ^^) i ^(a/ 32 ) ^2 ri/( a ^ 2 ) k( a /^) __ 

P nn ^2(002) ’ P 90 ^ 2 ( 002 ) ^^(0A0) £2 z0(OA 


1 , 0 ) 


to) 


(153) 

(154) 

(155) 

(156) 

(157) 

(158) 


(159) 


for j = In addition we obtain for j = z 

( a/ 9 l ) ^1 Tjrr(<*/ 31 ) n (a/ 32 ) ^2 tt/( q /^ 2 ) i jL( a £ 2 ) 

P Y^^ 001 ) ^ 12^3(001) - ^(0,0,1) + *«(0 A 1> 

’ * f v 1 

These relations provide 12 additional equations in the 2nd derivatives Whj mn y 

In summary, the imposition of traction continuity between neighboring cells and between the 
subcells of the cell in the radial, angular and axial directions provides a total of 44 equations in the 
2 nd time derivatives 
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4.13* Imposition of the Displacement Continuity Conditions. 

Continuity of displacements in the radial direction. The continuity conditions, imposed 
in the average sense, on the displacement vector at the inner surfaces as well as between neighboring 
cells in the r-direction, eqns (43) and (46), yield 


y {\ pl ) ^iu/C^T) , ^lry( 1 ^ 7 ) _ u/( 2 / 37 ) , y( 2 /? 7 ) _ fftpr/C 2 /^)] (p, 9 ,«) _ Q 

* 1(000) + 2 ”1(100) + 4 ”l(200) ” 1(000) + 2 ^1(100) 4 ”l(200)J ~ u 


(160) 


[w(W 4. , f|w( 2 ^ 7 )i(p,g,s) _ [*(1W _ , ^wW](p+l,?,s) 

i”l(000) + 2 ”^1(100) + 4 ”l(200)J ~ l” 1(000) 2 Wl (“») + 4 ” 1(200) J 

Iw( 1 ^ 7 ) . ^w( 107 ) 1 _ txA 2 /^) 1 ^ixt-I 2 /^) _ ^iy(W|(p,?,s) _ 

l 2(000) ^ 2 2 (100) + 4 **2(200) ^2(000) + 2 KK 2(100) 4 2(200)1 ~ 

[w(W , < kw&Pi) , ^ta^iOw,.) - “ ^ ^ 

1^2(000) + o ”2(100) + a ” 2(200) J 


0 


(161) 

(162) 

(163) 


Obviously, eqns (161) and (163), are not applicable at the boundary cell p = N p . These relations 
are replaced by the boundary conditions that are prescribed at the surface r = R\ of this cell as 
will be discussed in Section 4.14. 


Continuity of displacements in the angular direction. The continuity conditions, im- 
posed in the average sense, on the displacement vector, at the inner surfaces as well as between 
neighboring cells in the 0-direction, eqns (44) and (47), yield 
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(164) 
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Equations (165) and (167) are not applicable at the boundary cell q = N q . These relations are 
replaced by the boundary conditions that are prescribed at the surface 6 = © of this cell as will be 
discussed in Section 4.14. 


Continuity of displacements in the axial direction. The continuity conditions, imposed 
in the average sense, on the displacement vector in the periodic z-direction, eqns (45), provide 
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where e zz is the global axial strain that will be discussed in Section 4.15. This last relation follows 
from a homogenization procedure in the periodic direction that has been presented in Ref. [19]. 

These displacement continuity conditions provide altogether 44 relations in the 2nd time deriv- 
atives 

4 . 14 . Imposition of the Boundary Conditions. The tractions at the bottom (r = Rq) and 
top (r = i?i) surfaces are imposed by employing the relations: 



riAP i) i(l,g,s) jpO-Py) |(l,4,s) o_(l Pi) i(l.0. s ) 

^rj 1 r rj 1 u rj l f(0=-d( 1) /2 

(171) 


MWy) 1 (N p ,q,s) i p(W |0Vp,?,a)_ | (iV P’ 9 ’ s) 

r l 1 " r C? 1 rj l f (2) =<i ( W P ) /2 

(172) 

where j 
imposed 

= r y 9. By averaging these relations over the surfaces upon which these tractions are 
we obtain, respectively, 


ijAPy) _ jO-Pi) i(i,q,s) _ f (a j) 

(173) 


rr(2/?7) , d^Jip'y) i(N p ,q,s) _ / ■. /£ t \ 
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(174) 


(aQ'y) 

It follows, by employing the established expressions for the interfacial traction integrals m n ^ 
that the following equations are obtained in terms of the unknown 2nd time derivatives of the 
coefficients wj^y. 
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These relations readily imply the following equations: 
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(182) 
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If, on the other hand, time-dependent displacement is imposed at the boundary, the same 
averaging procedure is employed to establish the required expressions. For example, if the radial 
displacement u T is imposed at the right boundary 0 = 0, the following expression is obtained: 


K2, + yW®) + f = ■M’V) 


(185) 


4.15. Imposition of Plane Strain or Generalized Plane Strain Conditions . So far, we 
have established 104 relations in every cell for the determination of the 112 unknowns in 

this cell. The final set of 8 relations is determined from the imposition of either a plane strain or 
generalized plane strain condition in the periodic direction. 

By applying a homogenization procedure in the periodic z-direction, it can be shown that the 


following relation holds 


OW^ 1 } 

— 3(000) 

dz 


(186) 


for all a,/?, 7 = 1,2 in all the cells, where e Z2 is the far field average normal strain in the z- 
direction. Utilizing this relation, we can reduce the number of unknowns in each cell from 112 
to 104 by replacing the 8 unknowns with the single new unknown l zz . Consequently, the 

number of unknowns in each cell is 104, and the number of unknowns in the entire PGM composite 
comprised of N p N q cells is also 104A r p A r <J + 1. 

The global axial strain is related to the local strain = kF^oi) 35 f°fl° ws 

= E <% )h p% 4f T) (i87) 

P=1 a,/3 t 7=l 

where V = D(Rq + Ri)0(li + Z 2 )/2 is the total volume of the composite. 

Under plane strain condition this far field axial strain vanishes, namely, 

£** = 0 (188) 

For a generalized plane strain situation, on the other hand, the average normal stress in the 
z-direction vanishes: 

^.4ee t = o ( i8! » 

P=1 g=l a, 0,7= 1 

and the 2nd time derivative of this equation forms the required additional relation. 
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4.16. Summary of the Governing Equations. In summary, we have altogether 10AN p N q + l 
equations for the determination of the 2nd time derivatives of the l04N p N q 4- 1 displacement 
coefficients, , and the 2nd time derivative of the unknown global axial strain c zz . This 

system of equations can be represented in the following compact form 

AU(t) = f(i) + g(«) (190) 

where the structural stiffness matrix A contains information on the geometry and thermomechanical 
properties of the individual subcells (a/fy) within the cells comprising the functionally graded 
material. The displacement coefficient vector U(i) contains the 104 N p N q + 1 2nd time derivatives 
of the unknowns: 

U(t) = (()] (191) 

where 

OST’W = i.m,n = 0, 1,2. (192) 

The mechanical force vector f(t) contains information on the mechanical boundary conditions 
and the thermal loading effects generated by the applied temperature. In addition, the inelastic 
force vector g(t) appearing on the right-hand side of eqn (190) contains inelastic effects given in 
terms of the integrals of inelastic strain distributions that are represented by ^ n p 

The field variables that are expressed by the vector U (t) can be determined by integrating the 
above equation explicitly in a step-by-step timewise manner. To this end, let us introduce a time 
increment At. The time integration at time t yields 

U (t + At) = (At) 2 A -1 [f(t) + g(t)] + 2U(t) - U(t - At) (193) 

This difference expression which approximates the above 2nd order differential equation provides 
the displacement coefficients at the next time step t + At from the already known quantities at the 
current time t and the previous time step t — At (assuming that at time t = 0 the thermoelastic 
field distribution in the composite is known). Note that for time-independent material properties 
the matrix A has to be inverted just one time. This procedure is continued until the desired final 
time is reached. 


5. APPLICATION 

To demonstrate the potential of the outlined theoretical framework of the two-dimensional higher- 
order theory for cylindrical functionally graded materials with dynamic loading capability, a small 
computer code was developed to simulate one-dimensional wave propagation due to impulse loading 
in a layered half-plane as a special case. The half-plane consists of alternating layers of steel and 
polymeric material (PMMA) 0.025 and 0.0784 cm thick, respectively. Both materials are assumed 
to be linearly elastic with the relevant material parameters given in Table 1. 



Young’s modulus, E (dyne/cm 2 ) 

Mass density, p (gm/cm 3 ) 

Stainless steel 

1.258 x lO 17 

7.9 

PMMA 

0.089 x 10 12 

1.15 


Table 1. Material constants of a half-plane consisting of alternating steel and PMMA layers. 
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For an impacted linearly elastic composite by a spatially uniform loading applied normal to the 
layering, an exact solution by ray theory can be constructed, Ref. [23]. Figures 2 and 3 illustrate 
the predictions by the ray theory taken from this reference for the case when the half-plane is 
impacted by a unit impulsive normal stress applied at the center of the first layer (steel) at time 
t = 0 sec. The normal stress response with time is detected at the centers of the 11th, Fig. 2, 
and 12th layers, Fig. 3. The corresponding results generated by the computer code developed for 
this particular one-dimensional wave propagation case are included in the bottom portions of these 
figures. Excellent agreement between the exact ray-theory solution and the higher-order theory 
predictions is observed. 

6. PLANS FOR FUTURE WORK 

The completed generalization of the higher-order theory for cylindrical functionally graded materials 
with dynamic impact loading capability completely fulfills the objectives of the Grant NAG3-2411. 
As demonstrated in the above section, full benefit of the developed theoretical framework will be 
realized upon development of a general computer code enabling simulation of two-dimensional wave 
propagation in bi-directionally graded cylindrical structural components in the r — 6 plane. 
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8 . APPENDIX 


The nonzero strain coefficients in the Legendre polynomial expansion of the strain field 

subcell (a/3 7 ) of cell (p, q, s ) are given in terms of the displacement field microvariables by 
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The nonzero thermal stress coefficients in the Legendre polynomials expansion of the 

thermal stress in subcell (a 0 7) of cell (p, q, s ) are given in terms of of the temperature field mi- 
crovariables 

T(a0 7 ) _ r (aPy)r F (apy) 
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with similar expressions for and 
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Figure 1. A geometric model for the higher-order theory for cylindrical functionally graded mate 
rials/structures (H0TCFGM-2D). 
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Time (fis) 


Figure 2. Stress at the center of the 11th layer due to a unit step stress applied at the center of 
the first (top) layer: exact ray theory predictions from Ref. [23] (top figure); higher-order theory 
predcitions (bottom figure). The subscript 1 in <j\ indicates steel layer (layer 1) in the repeating 
unit cell sequence of the layered half-plane. 
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Time {\is) 


Figure 3. Stress at the center of the 12th layer due to a unit step stress applied at the center of 
the first (top) layer: exact ray theory predictions from Ref. [23] (top figure); higher-order theory 
predcitions (bottom figure). The subscript 1 in a 2 indicates PMMA layer (layer 2) in the repeating 
unit cell sequence of the layered half-plane. 
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